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Abstract 

The large deviation properties of equilibrium (reversible) lattice gases are mathematically 
reasonably well understood. Much less is known in non-equilibrium, namely for non 
reversible systems. In this paper we consider a simple example of a non-equilibrium 
situation, the symmetric simple exclusion process in which we let the system exchange 
particles with the boundaries at two different rates. We prove a dynamical large deviation 
principle for the empirical density which describes the probability of fluctuations from the 
solutions of the hydrodynamic equation. The so called quasi potential, which measures the 
cost of a fluctuation from the stationary state, is then defined by a variational problem for 
the dynamical large deviation rate function. By characterizing the optimal path, we prove 
that the quasi potential can also be obtained from a static variational problem introduced 
by Derrida, Lebowitz, and Speer. 
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1. Introduction 

In previous papers jHlU we have started the study of the macroscopic properties 
of stochastic non equihbrium systems. Typical examples are stochastic lattice gases 
which exchange particles with different reservoirs at the boundary. In these systems 
there is a flow of matter through the system and the dynamics is not reversible. 
The main difference with respect to equilibrium (reversible) states is the following: 
in equilibrium the invariant measure, which determines the thermodynamic proper- 
ties, is given for free by the Gibbs distribution specified by the Hamiltonian. On the 
contrary, in non equilibrium states the construction of the appropriate ensemble, 
that is the invariant measure, requires the solution of a dynamical problem. 

For equilibrium states, the thermodynamic entropy S is identified |Bl l20ll^ with 
the large deviation rate function for the invariant measure. The rigorous study 
of large deviations has been extended to hydrodynamic evolutions of stochastic 
interacting particle systems ^3 EI- Developing the methods of [Ci this theory has 
been extended to nonlinear hydrodynamic regimes jl5| . In a dynamical setting one 
may ask new questions, for example what is the most probable trajectory followed 
by the system in the spontaneous emergence of a fluctuation or in its relaxation 
to equilibrium. In the physical literature, the Onsager-Machlup theory [23 gives 
the following answer under the assumption of time reversibility. In the situation 
of a linear macroscopic equation, that is, close to equilibrium, the most probable 
emergence and relaxation trajectories are one the time reversal of the other. 

In [21 E] we have heuristically shown how this theory has to be modified for non 
equilibrium systems. At thermodynamic level, we do not need all the information 
carried by the invariant measure, but only its rate function S. This can be ob- 
tained, by solving a variational problem, from the dynamical rate function which 
describes the probability of fluctuations from the hydrodynamic behavior. The 
physical content of the variational problem is the following. Let p be the relevant 
thermodynamic variable, for instance the local density, whose stationary value is 
given by some function p{u). The entropy S{p) associated to some profile p{u) is 
then obtained by minimizing the dynamical rate function over all possible paths 
7r(i) — n{t,u) connecting p to p. We have shown that the optimal path TT*{t) is 
such that 7r*(—t) is a solution of the hydrodynamic equation associated to the time 
reversed microscopic dynamics, which we call adjoint hydrodynamics. This rela- 
tionship is the extension of the Onsager-Machlup theory to non reversible systems. 
Moreover, we have also shown that S solves an infinite dimensional Hamilton- J acobi 
equation and how the adjoint hydrodynamics can be obtained once S is known. 

In the present paper we study rigorously the symmetric one dimensional exclu- 
sion process. In this model there is at most one particle for each site of the lattice 
{—N,..., N} which can move to a neighboring site only if this is empty, with rate 
1/2 for each side. Moreover, a particle at the boundary may leave the system at 
rate 1/2 or enter at rate 7-/2, respectively 7+/2, at the site —N, respectively 
+N. In this situation there is a unique invariant measure p^ which reduces to a 
Bernoulli measure if 7_ =7+. On the other hand, if 7_ ^7+, the measure p^ 
exhibits long range correlations [7| [^ and it is not explicitly known. By using a 
matrix representation and combinatorial techniques, Derrida, Lebowitz, and Speer 
[HI [HI have recently shown that the rate function for p^ can be obtained solving 
a non linear boundary value problem on the interval [—1, 1]. We here analyze the 
macroscopic dynamical behavior of this system. The hydrodynamic limit for the 
empirical density has been proven in ^1 . We prove the associated dynamical 
large deviation principle which describes the probability of fluctuations from the 
solutions of the hydrodynamic equation. We then define the quasi potential via the 
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variational problem mentioned above. By characterizing the optimal path we prove 
that the quasi potential can also be obtained from a static variational problem 
introduced in |H1 E] . Using the identification of the quasi potential with the rate 
function for the invariant measure proven in we finally obtain an independent 
derivation of the expression for the thermodynamic entropy found in |8| |9j . 

2. Notation and results 

For an integer iV > 1, let An := [-N,N] nZ = {-N,...,N}. The sites of 
An are denoted by x, y, and z while the macroscopic space variable (points in 
the interval [—1,1]) by u, v, and w. We introduce the microscopic state space 
as Sat := {0,1}^" which is endowed with the discrete topology; elements of Sat, 
called configurations, are denoted by rj. In this way ri{x) S {0, 1} stands for the 
number of particles at site x for the configuration 77. 

The one dimensional boundary driven simple exclusion process is the Markov 
process on the state space Eat with infinitesimal generator 

Ln := i-,Ar + Lq^n + L+^N 

defined by 

x=~N 

{L±,Nf)iv) —h± + {l-l±H±N)][f{<J^^v)^.f{v)] 

for every function / : Sat — > K- In this formula a^'^r] is the configuration obtained 
from r] by exchanging the occupation variables r]{x) and 77(2/): 

(viy) a z^x 
ri{x) if z = y 
ri{z) li z^x,y 

and (T^r; is the configuration obtained from rj by flipping the configuration at x: 
{0-11) {z) := ?7(z)[l~4J + S,,S~v{z)], 

where 6x,y is the Kronecker delta. Finally 7± G (0,oo) are the activities of the 
reservoirs at the boundary of Aat. 

Notice that the generators are speeded up by N'^ ; this corresponds to the diffusive 
scaling. We denote by rjt the Markov process on J^n with generator L n and by P,, its 
distribution if the initial configuration is rj. Note that is a probability measure on 
the path space £'(R+, S^v), which we consider endowed with the Skorohod topology 
and the corresponding Borel cr-algebra. Expectation with respect to P^ is denoted 
byE^. 

Our first main result is the dynamical large deviation principle for the measure 
P,,. We denote by (•, •) the inner product in L2([— 1, 1], duj and let 

M -.^ {p e Loo{[-l,l],du) : 0<p(m)<1 a.e.} (2.1) 

which we equip with the topology induced by weak convergence, namely pn —>■ p in 
Ai if and only if {pn, G) {p, G) for each continuous function G : [— 1, 1] — > R; we 
consider M. also endowed with the corresponding Borel cr-algebra. Let us define 
the map tt^ : E Af ^ as 

^^(r?) := y ii{x)l{\- + — (2.2) 
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where 1{A} stands for the indicator function of the set A; namely (77) is 

the empirical density obtained from the configuration rj. Notice that tt^ {rj) G A4, 
i.e. < 7r^(?7) < 1, because ri{x) G {0, 1}. 

Let 77^ be a sequence of configurations for which the empirical density 77^(77^) 
converges in as t 00, to some function p, namely for each G G C([— 1, 1]) 

/■l^(t+2Tv) /■! 

lim (^^(77^),G) = lim ^ r;^(x) / duG{u)^ du p{u)G{u) 

(2.3) 

where we used the notation a A 6 := min{a, h] and aVb :— max{a, b}. If (|2.3|l holds 
we say that the sequence {r]^ : > 1} is associated to the profile p £ M. 

For T > and positive integers m, n we denote by G^'"([0,r] x [—1, 1]) the 
space of functions G: [0,T] x [—1, 1] ^ K. with m continuous derivatives in time, 
n continuous derivatives in space and which vanish at the boundary: G(-, ±1) = 0. 
Let also D{[0,T],M) be the Skorohod space of paths from [0,r] to M equipped 
with its Borel tr-algebra. Elements of D(\0,T], A4) will be denoted by n{t) = 
7r(t, u). 

Let p± := 7±/[l + 7±] G (0, 1) be the density at the boundary of [—1, 1] and fix 
a function p G which corresponds to the initial profile. For H G Gq'^([0,T] x 
[-1, 1]), let Jt.h,p = Jh ■■ D{\0, T],M) — > M be the functional given by 

JniT^) {n{T),H{T))-{p,H{Q))- £dt(TT{t),dtH{t) + ^AH{t)) 

-l£dt{xinit)),{VHit))'), (2.4) 

where V denotes the derivative with respect to the macroscopic space variable 
u, A is the Laplacian on (—1, 1), and we have set x(^) '^(1 ~ '^)- Let finally 
/t( • \p) ■■ D{[0, T],M) — y [0, +00] be the functional defined by 

/t(^Ip) := sup Jh{tt) . (2.5) 

HeCo^''([O.T]x[-l,l]) 

Notice that, if n(t) solves the heat equation with boundary condition 7r(i, ±1) = p± 
and initial datum 7r(0) = p, then It{t^\p) = 0. 

Theorem 2.1. Fix T > and a profile p (z A4 bounded away from and 1, namely 
such that there exists S > with S < p < 1 — S a.e. Consider a sequence rj^ of 
configurations associated to p. Then the measure P^n o (tt^)^^ on D{\0,T], Ai^ 
satisfies a large deviation principle with speed N and convex lower semi continuous 
rate function It{-\p)- Namely, for each closed set C C D{[0,T], M) and each open 
set O C D{[0,T],M), 

limsup 4 logP„" [tt^ G C] < - inf It{tt\p) 

N^oo N Tree 

hminf llogP,«[7r^ G O] > - inf /T(7r|p) . 

It is possible to obtain a more explicit representation of the functional It{-\p), 
see Lemma 13.61 below. If the particle system is considered with periodic boundary 
conditions, i.e. An is replaced by the discrete torus of length N, this Theorem has 
been proven in jl7|. As we shall see later, the main difference with respect to the 
case with periodic boundary condition is the lack of translation invariance and the 
fact that the path 7r(t, •) is fixed at the boundary. 
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We now define precisely the variational problem mentioned in the introduction. 
Let p & M he the hnear profile p{u) :— [p_ (1 — u) + p+(l + m)] /2, u G [—1,1], 
which is the density profile associate to the invariant measure /i^ , see Section 01 
below. We then define V : Ai ^ [0, +oo] as the quasi potential for the rate function 

V(p) := inf inf IrMp) . (2.6) 

T>0 7r(.) : Tr{T)=p 

which measures the minimal cost to produce the profile p starting from p. 

Let us first describe how the variational problem H2.6|l is solved when 7_ = 7+ = 
7. In this case p = 7/(1 + 7) is constant and the process is reversible with respect 
to the Bernoulli measure with density p. We have that It{t^\po) = if 7r(i) solves 
the hydrodynamic equation which for this system is given by the heat equation: 

dtpit) - (l/2)Ap(t) , 

p(t,±l) = p±, (2.7) 
p(0,-)=Po(-) . 

Note that p{t) ^ p as t ^ 00. 

It can be easily shown that the minimizer for (|2.6|) , defined on the time interval 
(—00, 0] instead of [0, +00) as in H2.6|l . is given by Tr*{t) = p{—t), where p{t) is the 
solution of H2.7|l with initial condition po — p. This symmetry of the relaxation and 
fluctuation trajectories is the Onsager-Machlup principle mentioned before. 

Moreover the quasi potential V{p) coincides with the entropy of the Bernoulli 
measure with density p, that is, understanding OlogO = 0, 

V{p) = Soip) := du \piu) log ^ + [1 - piu)] log (2.8) 
7-1 L p 1-p J 

In the context of Freidlin-Wentzell theory ^1] for diffusions in M" , the situation 
just described is analogous to the so called gradient case in which the quasi po- 
tential coincides with the potential. This structure reflects the reversibility of the 
underlying process. In general for non gradient systems, the solution of the dynam- 
ical variational problem, or of the associated Hamilton-Jacobi equation, cannot be 
explicitly calculated. The case 7+ ^ 7_ is analogous to a non gradient system, but 
for this particular model we shall prove that the quasi potential V{p), as defined in 
()2.6|l . coincides with the functional S{p) defined by a time independent variational 
problem introduced in [HI E] which is stated below. This is the second main result 
of this paper. 

Denote by C^([— 1,1]) the space of once continuously differentiable functions 
/ : [—1, 1] R endowed with the norm H/Uci ■= sup,^gj_]^ {l/(''^)l + l/'(")l}- Let 



{/eCi([-l,l]) : /(±l) = p±, [p+-p^]riu)>0,ue[-l,l]] , (2.9) 

where /' denotes the derivative of /. Note that if / e then < p_ A /9+ < f{u) < 
P-V p+ < 1 for all -1 < M < 1. 
For p £ M and / e we set 

GipJ):^ /^4p(.)logf^+[i_,(„)]logi:iZH+log. 
J-i L f{u) l-/(w) 



[p+-p-]/2 



(2.10) 



and 



S{p) := sup Gip J) . (2.11) 

Theorem l4 . SI below . which formalizes the arguments in 0, states that the supremum 
in H2.11() is uniquely attained for a function / which solves a non linear boundary 
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value problem. We shall denote it by F = F{p) to emphasize its dependence on p; 
therefore Sip) = G {p, F{p)) . 

Theorem 2.2. Let V and S as defined in \2. 6]) and K2.11^ . Then for each p E Ai 
we have V{p) = S{p). 

In the proof of the above theorem we shall construct a particular path n*{t) in 
which the infimum in H2.6|l is almost attained. As recalled in the introduction, by 
the heuristic arguments in 0j, Tr*{—t) is the solution of the hydrodynamic equation 
corresponding to the process with generator L^, the adjoint of Xat in L2{'^Ni dp,^) 
and initial condition p. In analogy to the Freidlin-Wentzell theory jj^, we expect 
that the exit path from a neighborhood p to a neighborhood of p should, with 
probability converging to one as iV "f oo, take place in a small tube around the path 

The optimal path can be described in a rather simple fashion. Recalling that 
we denoted by _F = F{p) the maximizer for H2.11() . consider the heat equation in 
[— 1, 1] with boundary conditions p± and initial datum F: 

dMt) = (l/2)A<i>(0 , 

$(t,±l)-p±, (2.12) 
$(0, •) - . 

We next define p*{t) — p*{t,u) by 

p*{t) m + ^t)[l-m]-^^ ■ (2.13) 

In view of H4.3|l below, p*(0) = p and, by Lemma 15.61 limt^oo P*{t) = P- The 
optimal path 7r*(t), defined on the time interval (— oo,0] instead of [0, +oo) as in 
(ITHji . is then given by 7r*(t) = p*{-t). 

From the dynamical large deviation principle we can obtain, by means of the 
quasi potential, the large deviation principle for the empirical density when the 
particles are distributed according to the invariant measure of the process r]t ■ Note 
that the finite state Markov process -qt with generator L^r is irreducible, therefore 
it has a unique invariant measure pJ^ . 

Let us introduce Vn '■= ° i''^^) which is a probability on M and describes 
the behavior of the empirical density under the invariant measure. In 7. 12, 13il24| 
it is shown, see also Section O below, that Vn satisfies the law of large numbers 
Vn =^ 5p in which =4> stands for weak convergence of measures on M. and p is the 
linear profile already introduced. 

Since p is globally attractive for (|2.7(l . the quasi potential with respect to p de- 
fined in (|2.6|) gives the rate function for the family Vn- In 0I5 we have heuristically 
derived this identification via a time reversal argument. For the present model a 
rigorous proof, in the same spirit of the Freidlin-Wentzell theory, is given in 0; 
that is we have the following theorem. 

Theorem 2.3. Let V as defined in \2.b]) . Then the measure Vn satisfies a large 
deviation principle with speed N and rate Junction V . 

The identification of the rate function for Vn with the functional S now follows 
from Theorems O O and lO 



Corollary 2.4. Let S as defined in \2.11]) . The measure Vn satisfies a large 
deviation principle on M with speed N and convex lower semi continuous rate 
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function S . Namely for each closed set C d A4 and each open set O C A4, 
limsup 4 log M^k" eC]<- inf S{p) 
liminf 1 log fi^'in^ e O] > - inf Sip) 

As already mentioned, the rate function 5 has been first obtained in [SJ |5] by 
using a matrix representation of the invariant measure p^ and combinatorial tech- 
niques. By means of Theorems 12.11 12.21 and 12.31 we prove here, independently of 
[HI El J the large deviation principle by following the dynamical/ variational route ex- 
plained in 0] which is analogous to the Freidlin-Wentzell theory for diffusions 
on R". 

We remark that it should be possible, modulo technical problems, to extend 
Theorems 12.11 and 12.31 to other boundary driven diffusive lattice gases, see 0] for 
a heuristic discussion. The characterization of the rate function for the invariant 
measure as the quasi potential allows to obtain some information on it directly 
from the variational problem (|2.6|) . In particular, in Appendix fXl we discuss the 
symmetric simple exclusion in any dimension and get a lower bound on V in terms 
of the entropy 5*0 of the equilibrium system. In the one dimensional case, this 
bound has been proven in jH| EI by using instead the variational problem 1)2.11(1 . 

Outline. In Section |31 we recall the hydrodynamic behavior of the boundary driven 
exclusion process and prove the associated large deviation principle described by 
Theorem 12.11 In Section 0] and which are more technical, we state and prove 
some properties of the functional S which is then shown to coincide with the quasi 
potential V. Finally, in Appendix fXl we consider the symmetric simple exclusion 
in any dimension and prove a lower bound on V. 

3. Dynamical behavior 

In this section we study the dynamical properties of the empirical density for 
the boundary driven simple exclusion process in a fixed (macroscopic) time interval 
[0,T]. In particular we review the hydrodynamic limit (law of large numbers) and 
prove the corresponding large deviation principle. This problem was considered 
before by Kipnis, OUa and Varadhan in (17) for the exclusion process with periodic 
boundary condition. For this reason, we present only the modifications needed in 
the argument and refer to ^| ^1 12] for the missing arguments. 

As already stated, the invariant measure p^ is not known explicitly but some 
of its properties have been derived. For example, the one site marginals or the 
correlations can be computed explicitly. To compute the one site marginals, which 
will be used later, let p^ [x) = E^n \ri{x)\ for —TV < x < N. Since p^ is invariant, 
E^N[LNrj{x)] — 0. Computing Lpfijix), we obtain a closed difference equation for 
p^ix): 

' (AAr/9^)(x) = for -iV+ 1 < x < A^- 1 , 

p^{N - 1) - p^{N) + 7+[l - P^'iN)] - p^'iN) = , 
p^i^N + 1) - p^i-N) + 7- [1 - p^i-N)] - p^i-N) = . 

In this formula, Ajv stands for the discrete Laplacian so that (Ajv/)(a;) = f{x + 
1) -I- f{x — 1) — 2f{x). The unique solution of this discrete elliptic equation gives 
the one-site marginals of p^ . 

Denote by = the product measure on Ejv with marginals given by 

: rjix) = 1} = p^{x) 
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and observe that the generators L-,7v, i+,Ar are reversible with respect to . 

Denote by {t^ : x £ Z} the group of translations in {0, 1}^ so that {txC){z) = 
C(x + z) for all a;, z in Z and configuration ( in {0, 1}^. Translations are extended 
to functions and measures in a natural way. Eyink, Lebowitz and Spohn '12' and 
De Masi, Ferrari, laniro and Presutti proved that 

for every local function / and u in (—1,1). Here p is the unique solution of 

(l/2)Ap = 0, 
P(±l) = P± , 

namely p is the linear interpolation between p_ and p+ and {i^a:0<Q;<l} stands 
for the Bernoulli product measure in {0, 1}^ with density a and p± = 7±/[l + 7±] 
is the density at the boundary of [—1,1]. 

3.1. Hydrodynamic limit. Recall that, for each configuration r] G S^r, we denote 
by TT^ — TT^ {rf) G A4 the empirical density obtained from r], see equation (|2.2|l . We 
say that a sequence of configurations {rj^ : > 1} is associated to the profile 7 if 
(|2.3|) holds for all continuous functions G : [— 1, 1] ^ M. The following result is due 
to Eyink, Lebowitz and Spohn |13| . 

Theorem 3.1. Consider a sequence rj^ associated to some profile po £ A4. Then, 
for all t > Q, {t) = TT^ {rjt) converges (in the sense ) in probability to p{t, u), 
the unique weak solution of 

dtp = (l/2)Ap , 

p(i,±l)=p±, (3.1) 
p(0,-) = Po(-) . 

By a weak solution of the Dirichlet problem (|3.1|l in the time interval [0,T], we 
understand a bounded real function p which satisfies the following two conditions. 

(a): There exists a function A{t,u) in L^([— 1,1] x [0, T]) such that 
ds [ dup{s,u)(\7H){u) 



J-i 



{p+H{l)-p^H{-l)}t - ^ dsj^^ 



du A(s, u)H{u) 



for every smooth function H : [—1, 1] ^ M and every < t < T . A{t,u) 
will be denoted by {\/p){t,u). 
(b): For every function H : [—1, 1] ^ M of class C^([— 1, 1]) vanishing at the 
boundary and every < t < T, 

1 .1 
dup{t,u)H{u) — I dupo{u)H{u) 



-1 

{1/2) j^dsj du{\7p)[s,u)[\IH){u) . 



The classical estimates gives uniqueness of weak solutions of equation H3.1|l . 
Note that here the weak solution coincides with the semi-group solution p{t) = 
p + e*^ ^^(po — p)j where p is the stationary profile and A° is the Laplacian with 
zero boundary condition. 
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3.2. Super-exponential estimate. Wc now turn to the problem of large devia- 
tions from tlic liydrodynamic limit. It is well known that one of the main steps in 
the derivation of a large deviation principle for the empirical density is a super- 
exponential estimate which allows the replacement of local functions by functionals 
of the empirical density in the large deviations regime. Essentially, the problem 
consists in bounding expressions such as (V,/^)^™ in terms of the Dirichlet form 
(— Ljv/, /)^«. Here ^ is a local function and (•, indicates the inner product 
with respect to the invariant state /i^. 

In the context of boundary driven simple exclusion processes, the fact that the 
invariant state is not known explicitly introduces a technical difficulty. Following 
|19| we fix , the product measure defined in the beginning of this section, as 
reference measure and estimate everything with respect to . However, since ly^ 
is not an invariant state, there are no reasons for {—L^f, f)^N to be positive. The 
first statement shows that this expression is almost positive. 

For a function / : S^r R, let 

JV-l „ 
x=-N •' 

Lemma 3.2. There exists a finite constant Co depending only on ^± such that 

(io,JV < -—DNif) + CoNifJ),. 
for all functions f : E^r — s- R 

The proof of this lemma is elementary and left to the reader. Notice on the 
other hand that both (— L+^jv/, , {—L-.Nf:f)w^ a-re positive because is a 
reversible state by our choice of the profile . 

This lemma together with the computation presented in [21 p. 78] for non- 
reversible processes, permits to prove the super-exponential estimate. The state- 
ment of this result requires some notation. For a cylinder function VE", denote the 
expectation of ^ with respect to the Bernoulli product measure Va by \E'(a): 

For a positive integer i and —N<x<N, denote the empirical mean density on a 
box of size 2£ + 1 centered at x by ri^{x), namely 

^ Jam ^ '^(y) ' 

where Ai{x) = AN,i{x) ^ {y e An : \y - x\ < £}. Let H e C([0, T] x [-1, 1]) and 
a cylinder function. For £ > 0, define also 

V"f{t,v) = ^ 5] i?(t,x/7V){r.*(77)- ^(77^^(0:))}, 

X 

where the summation is carried over all x such that the support of t^"^ belongs to 
Aat. For a continuous function G : [0,T] R, let 



^ f dsGis)[r]si±N)-p±] 
Jo 



Theorem 3.3. Fix H in C{[0,T] x [-1, 1]), G G C([0,T]), a cylinder function 
and a sequence {rj^ G T,^ ■ N > 1} of configurations. For any S > we have 



lim sup lim sup — log P,^^ 



^ vS;^{t,ijt)dt 



> S 



= —00 , 



limsup — logP jv 



\W^\ > 6 



-oo 



3.3. Upper bound. The proof of the upper bound of the large deviation principle 
is essentially the same as in There is just a slight difference in the definition 
of the functionals Jh due to the boundary conditions. 

For H in Cg'^([0, T] x [—1, 1]) consider the exponential martingale defined 

by 



H 



expi N 



NJo 



-Ni.-is)M(s)) ^ ^2 ^^^^Ni.- is) Mis)) 



An elementary computation shows that 

M# = expiv{j^^(7r^*t,) + V^_, + ^^(e)} , 

where limg^o C'ffC^) = 0, stands for the approximation of the identity 
{2e)^^l{u G [— e,e]}, * stands for convolution, 



N,6 



VHi--l) 



and 4-0(7/) =7/(0)[l-?7(l)]. 

Fix a subset A of D{[0,T],M) and write 



1 



log I 



1 



■logE,- 



m^{m^)-H{7t'^ e A} 



N ° " ' ' N 

Maximizing over tt^ in A, we get from previous computation that the last term is 
bounded above by 



1 



CH{e) . 



- inf JH(7r*te) + — logE i^ 

Denote by the measure P,jjv Since the martingale is bounded by exp{CA^} 
for some finite constant depending only on H and T, Theorem 13.31 holds for 
in place of P^n . In particular, the second term of the previous formula is bounded 
above by Ch{£,N) such that limE^o limsup^^oQ Ch{£,N) — 0. Hence, for every 
e > 0, and every H in Co'^([0,T] x [-1, 1]), 

lim sup ^ log P^N [tt^ e A] < - inf JniTr * te) + Cfjie) , 



N 



TreA 



where limg^o — 0. 

Assume now that the set A is a compact set K. Since Jh{- * i-e) is continuous 
for every H and e > 0, we may apply the arguments presented in Lemma 11.3 of 
|25| to exchange the supremum with the infimum. In this way we obtain that the 
last expression is bounded above by 

limsup4logIP»,«[7i'^ K] < - inf sup | JH(7r * i^) + C^(e)| . 

Letting first £ i 0, since Jh{t^ * t^) converges to J//(7r) for every H in Cq'^([0, T] x 
[—1,1]), in view of the definition (|2.5|l of /T(7r|7), we deduce that 

limsup^logP iv[7r^ e AT] < - inf /T(7r|7) , 

which proves the upper bound for compact subsets. 
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To pass from compact sets to closed sets, we have to obtain "exponential tight- 
ness" for the sequence P^jv [tt^ S •] . The proof presented in Jj, for the non inter- 
acting zero-range process is easily adapted to our context. 

3.4. Hydrodynamic limit of weakly asymmetric exclusions. Fix a function 
H in Cq'^([0; T] X [—1, 1]) and recall the definition of the martingale . Denote by 

the probabihty measure on ^([O, T], Eat) defined by P,^„ [A] = E,,jv [M^1{A}]. 
Under P^, the coordinates {rjt < t < T} form a Markov process with generator 

= L+^N + L^^M + L-,N, where 

{L^.NfM - — E e-{^(*'[-+il/^)-^(*'-/^)H''(-+i)-''(-)>[/(a-'-+i77)-/(7y)] . 

The next result is due to Eyink, Lebowitz and Spohn J^l- Recall x{p) = p(l^p)- 

Lemma 3.4. Consider a sequence rj^ associated to some profile j (E Ai and fix H 
in Cg'^([0, T] X [—1, 1]). Then, for all t > 0, 7r^(i) — {rjt) converges in probability 
(in the sense ^2.S\) ) to p{t,u), the unique weak solution of 

9tP=(l/2)Ap~V{x(p)Vi7}, 

p{t,±l)=p±, (3.2) 
p(0,-)=7(-)- 

As in subsection 13.11 by a weak solution of the Dirichlct problem (|3.2|) in the 
time interval [0,T], we understand a bounded real function p which satisfies the 
following two conditions. 

(a) : There exists a function A{t,u) in L^{[—1, 1] x [0,r]) such that 

ds J dup{s,u){S/G){u) 

= {p+G{l) - p-G{-l)}t - J^dsJ duA{s,u)G{u) 

for every smooth function G: [—1, 1] ^ M and every < t <T. A{t, u) will 
be denoted by {S/p){t,u). 

(b) : For every function G E G^ ([^1' 1]) vanishing at the boundary and every 
t > 0, 

"1 

du p{t,u)G{u) — / duj{u)G{u) — 
1 J-i 

ds j ^du {VG){u)[ - (l/2)(Vp)(5, u) -h x(p(s, u)){VH){s, u)] . 

The classical H^i estimates gives uniqueness of weak solutions of equation H3.2|l . 

3.5. The rate function. We prove in this subsection some properties of the rate 
function It{ ■ I7). We first claim that this rate function is convex and lower semi 
continuous. In view of the definition of It{ ■ I7), to prove this assertion, it is enough 
to show that Jh is convex and lower semi continuous for each H in Cg'^([0,T] x 
[—1, 1]). It is convex because x(a) = a(l — a) is a concave function. It is lower semi 
continuous because for any positive, continuous function G : [0,r] x [—1, 1] ^ M 
and for any sequence tt" converging to tt in D{[Q,T],Ai), 

dt (x(7r(t)), G(0) - lim r dt {x{n{t) * i,), G{t)) 

= lim lim / (x(7r"(<) * te),G(<)) . 
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Since % is concave and G positive, a change of variables shows that this expression 
is bounded below by 

limlimsup / dt {x{n''{t)),G{t) * L,) = limsup [ dt {x{'K''{t)),G{t)) 

because G is continuous and x is bounded. This proves that Jr is lower semi 
continuous for every H in Cq'^([0,T] x [—1, 1]). 

Denote by the subset of D{[0,T\,M.) of all paths TT{t,u) which satisfy the 
boundary conditions 7r(0, •) = 7(-), 7r(-,±l) = p±, in the sense that for every 
< to < ii < T, 

j-ti r — l+S 

lim / dt - Tr{t, u) du = {ti — to)p- 

and a similar identity at the other boundary. 

Lemma 3.5. It {ttIj) = oo if n does not belong to D^. 

Proof: Fix tt in D{[0, T],M) such that Iri-^ll) < oo. Wc first show that 7r(0, •) = 
7(-). For 6 > 0, consider the function Hs{t,u) = hs{t)g{u), hs{t) = (1 — 5~^t)^, 
g{-) vanishing at the boundary ±1. Here o+ stands for the positive part of a. Of 
course, Hs can be approximated by smooth functions. Since n is bounded and since 
t Tr(t, •) is right continuous for the weak topology, 

lim-^^ff^W = (7r(0),ff) - (7,5) , 

which proves that 7r(0) = 7 a.s. because /T(7r|7) < 00. 

A similar argument shows that 7r(t, ±1) = p±; to prove this statement we may 
consider the sequence of functions Hs{t,u) = h{t)gs{u), where h{t) approximates 
the indicator of some time interval [toi^i] and where 



A- {A + b){l + u)/6 if-1 < u < -1 + 5, 
-b if-l + 5<u<l. 



Here A > is large and fixed and b = b{A, (5) > is chosen for the integral over 
[— 1, 1] of g'g to vanish. □ 

Fix TT in D-y and denote by Hi(7r) the Hilbert space induced by Cq^{[0,T] x 
[— 1, 1]) endowed with the inner product (•, •)7r defined by 

{H,G)„ = £ dt j'^dux{n){VG){VH) . 

Lemma 3.6. Fix a trajectory tt in and assume that /T(7r|7) is finite. There 
exists a function H in Hi (tt) such that w is the unique weak solution of 

dt-K = (l/2)A7r - V{x(7r)7r(l - ■k)VH} , 
7r(t,±l)=p±, (3.3) 
7r(0,-)=7(-)- 



Moreover, 



/T(7r|7) = i^/^)£dt j'^dux{7T){WHf . (3.4) 
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Wc refer the reader to El for a proof. One of the consequences of this 
lemma is that every trajectory t t—f 7r(i) with finite rate function is continuous in 
the weak topology, tt S C([0,T];7W). Indeed, by the previous lemma, for tt such 
that /T(7r|7) < oo, and every G in Cq([— 1, 1]), 

(7r(0,G)-(^(s),G) = (1/2) y*dr(7r(r), AG) +y'dr(x(^(r)),VGViJ) 

-(1/2){(VG)(1) p+ - (VG)(-1) p_}(t - s) 

for some H in Hi (tt) . Since G is smooth and H belongs to Hi (vr) , the right hand 
side vanishes as |t — s| ^ 0. 

3.6. Lower bound. Denote by the set of trajectories tt in D{[0,T], Ai) for 
which there exists H in Go'^([0,T] x [—1, 1]) such that it is the solution of H3.3|l . 
For each tt in D!^, and for each neighborhood A/"^ of tt 

liminf i-logP„«[7r~ e AA,] > -/t(^|7) ■ 

This statement is proved as in the periodic boundary case, see ^Hl- To com- 
plete the proof of the lower bound, it remains to show that for every trajectory tt 
such that /T(7r|7) < oo, there exists a sequence tt^ in such that lim/j ttj; = tt, 
limfc/T(7rfc|7) = /T(7r|7). 

This is not too difficult in our context because the rate function is convex and 
lower semi continuous. We first show that any path tt with finite rate function 
can be approximated by a path which is bounded away from and 1. Fix a path 
TT such that /T(7r|7) < oo. Fix S > and denote by p{t,u) the solution of the 
hydrodynamic equation (|3.1|l with initial condition 7 instead of pQ. Let tts = Sp + 
(1 — S)tt. Of course, tts converges to tt as (5 | 0. By lower semi continuity, /T(7r|7) < 
liminfa^o ^T(7r5|7). On the other hand, since /t(- I7) is convex, IxiT^sll) < (1 — 
5)It{t^\j) because p is the solution of the hydrodynamic equation and /t(p|7) = 0. 
This shows that limg^oTTs = tt, limg^o It {tts\j) = I{t^)- Since < 7 < 1, < 
p± < 1, TT^ is bounded away from and 1, proving the claim. 

Fix now a path tt with finite rate function and bounded away from and 1. We 
claim that this trajectory may be approximated by a path in D^. Since It{tt\^) < 
00, by Lemma l3.6l there exists H in Hi{tt) satisfying 1)3. 3|1 . Since tt is bounded away 
from and 1, Hi{tt) coincides with the usual Sobolev space Hi associated to the 
Lebesgue measure. Consider a sequence of smooth functions : [0, T] x [— 1, 1] 
M vanishing at the boundary and such that V-ff„ converges in L^([0,T] x [—1,1]) 
to VH . Denote by tt" the solution of 1)3.2(1 with instead of H . We claim that 
lim^^oo tt" = TT, lim„^oo /T(7r"|7) = /T(7r|7)- 

The proof that tt" converges to tt is divided in two pieces. We first show that 
the sequence is tight in G([0,T],A^) and then we prove that all limit points are 
solution of equation (|3.2|l . We start with a preliminary estimate which will be 
needed repeatedly. Recall that p is the stationary profile. Computing the time 
derivative of ^^-^ du{TT"'{t) — pY , we obtain that 

dt J du(V7r"(t))2 < C (3.5) 

for some finite constant independent of n. 

From the previous bound and since 7r"(i, u) belongs to [0, 1], it is not difficult to 
show that the sequence tt" is tight in G([0,T], A^). To check uniqueness of limit 
points, consider any limit point /? in C{[0,T], A4). We claim that /? is a weak 
solution of the equation (|3.2|l . Of course /3 is positive and bounded above by 1. 
The existence of a function A{s, u) in L^([— 1, 1] x [0, T]) for which (a) holds follows 



13 



from l|3.5|l . which guarantees the existence of weak converging subsequences. The 
unique difficulty in the proof of identity (b) is to show that for any < t < T, G 



m 



L2([o,r] X [-1,1]), 



hm f ds{xi7r'\s)),G{s)) = f ds {x{Pis)),G{s)) (3.6) 



n — *oo 



for any sequence tt" converging to (3 in C{[0,T], A4) and satisfying 1)3.5(1 . This 
identity holds because for any 5 > 

lim f ds{xi7r^is)*Ls),G{s)) = f (x(/9(s) * i^), G(s)) 
"-'°°7o Jo 

and because, by Schwartz inequality and |x(a) ~ x{b)\ < |a — 

''*ds(x(^"(s)*.,-)-x(^"(s)),G(s))^' 

< /* ds (G(s)2) /* ds ([7r"(s) * is - ^"(s)]^) . 
Jo Jo 

It is not difhcult to show, using estimate (|3.5|l . that this term vanishes as (5 J, 0, 
uniformly in n, proving 1(3.6(1 . In conclusion, we proved that the sequence tt" is 
tight in C{[0,T], Ai) and that all its limit points are weak solutions of equation 
((3.2() . By uniqueness of weak solutions, tt" converges in C{[0,T],M) to n. 

It remains to see that It{t^"'\^) converges to /T(7r|7). Since tt" ^ tt and It{ • I7) 
is lower semi continuous, we just need to check that lim sup„ It (tt"] 7) < It{tt\j). 
Here again the concavity and the boundness of x help. Since converges in 

to V-ff and x is bounded, the main problem is to show that 



limsup / dt ixiTT- it)), {V H {t)y) < / dt{xi7r{t)),iVH{t))') . 

n^oc Jo Jo 

Since tt * converges almost surely to tt as (5 | 0, 

^ dt{x{7r{t)),{S/H{t))') = r dt{x{n{t)*ts)A^H{t))^) 



-^-ojo 



T 



= lim lim / dt{x{TT"'{t)*Ls)A'^H{t))^) . 
Since x is concave, the previous expression is bounded below by 
lim limsup / dt (x(7r"(t)) * is, {\7H{t))^) . 

S—^O n^oc Jo 

Since x is bounded and (ViJ)^ integrable, a change of variables shows that the 
previous expression is equal to 

limsup r dt{x{7r"{t)),{S/H{t))') , 
concluding the proof of the lower bound. 



4. The rate function for the invariant measure 

In this section we discuss some properties of the functional S{p) which are needed 
later. The results stated here are essentially contained in [Hj, but, for the sake of 
completeness, we review them and give more detailed proofs. Without any loss 
of generality, from now on we shall assume that < p_ < p+ < 1. Recall the 
definitions of the set JF, 1(2. 9() . and of the functional G{p,f), 1(2. 10() . 
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The Euler-Lagrange equation associated to the variational problem H2.11|l is 
given by the non linear boundary value problem 

F" =(p-F) ) ' , in (-1,1), , , 

F(±l) = p± . 

We introduce the notation, which we will use throughout this section, 

niu) ^ n{p,F-u) = {p{u) - F{u)) J'^""^ . (4.2) 

F{u){l~F{u)) 

Using this notation equation (|4.1|l takes the form 

F" = F'n in (-1,1), 



F(±l) = p± 



(4.3) 



In order to state and prove an existence and uniqueness result for F € we 
formulate (|4.3(l as the integro-differential equation 

/ dv exp < / dw TZ{p, F; w) > 

F{u) = p_ + (p+ - p^y-^ 1 . (4.4) 

/ dv exp < / dw TZ{p, F; w) > 

We will denote its solution by F = F{p) to emphasize its dependence on p. We 
observe that if p = p then F = F{p) — p solves H4.4II and (|4.3|l . 

Notice that if F e C^([— 1, 1]) is a solution of the boundary value problem (|4.3|l 
such that F'{u) > for u e [—1,1], then F is also a solution of the integro- 
differential equation H4.4() . Conversely, if F e C^([— 1,1]) is a solution of H4.4|l . 
then F'{u) > 0, F"{u) exists for almost every u and H4.3|l holds almost everywhere. 
Moreover, if p e C(\-\, 1]), then F £ C^([-l, 1]) and H4.3|l holds everywhere. 

Remark 4.1. There are non monotone solutions of equation i4.S{ ). For example, 
for the constant profile p = 1/2, it is easy to check that the functions 

Fiu) = \ [l + sin(Au + (^)] ; 

satisfy equation {4-!^ for countably many choices of the parameters X and ip (fixed in 
order to satisfy the boundary conditions in However only one such function 

is monotone. In fact, under the monotonicity assumption on F, we will prove 
uniqueness (and existence) of the solution of the boundary value problem {4.3)) . 

The following theorem gives us the existence and uniqueness result for H4.4|l to- 
gether with a continuous dependence of the solution on p. Recall that we denote by 
C^([— 1, 1]) the Banach space of continuously differentiable functions / : [—1, 1] 
R endowed with the norm ll/Hc-i ■= suPue[-i,i] {1/(^)1 + l/'(^)l}- 

Theorem 4.2. For each p £ Ai there exists in T a unique solution F = F{p) of 
\4-4^ - Moreover: 

(i) z/p € C([— 1, 1]), then F = F[p) G C^([— 1, 1]) and it is the unique solution 
inTC]C^{[-l,l]) of U^ ; 

(ii) if pn converges to p in M as n oo, then Fn = F{pn) converges to 
F = F{p) m Ci([-l,l]),- 

(iii) fix T > and consider a function p = p{t,u) G C^''^([0,T] x [—1,1]). Then 
F - Fit,u) = F{p{t, ■)){u) e Ci'2([0,r] X [-1, 1]). 
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The existence result in Theorem 14 . 21 will be proven by applying Schauder's fixed 
point theorem. For each p & M. consider the map /Cp : ^ ([— 1, 1]) given by 

/ dveyi'pl I dwTZ{p, f;w) > 

/Cp(/)(u) := p_ + (p+ - p^y-jl L ■ (4.5) 

/ dv cxp < / dwTZ{p, f;w) > 

Let us also define the following closed, convex subset of Ci([-l,l]): 

B:={fe C\[-l, 1]) : /(±1) ^p±,b< f{u) <B}cT, (4.6) 

where, recalling we are assuming 7_ < 7+, 

^ P+ - P- 1- Q .^ P+- P- 7+ 
27+' ■ 2 7_ ■ 

Lemma 4.3. For each p G Ai, ICp is a continuous map on T and K,p{!F) C B. 
Furthermore JCp(B) has compact closure in C^([— 1, 1]). Hence, by Schauder's fixed 
point theorem, for each p & M. equation {4-4^ has a solution F = K,p{F) S B. 
Moreover, there exist a constant C G (0,oo) depending on p± such that for any 
p € M and any u,v € [—1, 1] we have \F'{u) — F'(v)\ < C \u — v\. 

Proof: It is easy to check that ICp is continuous and /Cp(/)(±1) = p±. Let us define 
5p := ^p(/)> we have 



g'^{u) ^ {p+ - p^)— ^-^^^^ (4.7) 



Since p{w) — f{w) < 1 — f{w), p{w) — f{w) > —f{w), and /'(w) > 0, we get 

(1-/)' ^ /' 
.1^ il. < 7^ < 

1-/ - - / 

which implies b < g'p{u) < B for all u G [—1,1]- In particular lCp{!F) C B. 

To show that JCp{B) has a compact closure, by Ascoli-Arzela theorem, it is 
enough to prove that g'^ is Lipschitz uniformly for f € B. Indeed, by using H4.7|l . it 
is easy to check that there exists a constant C — C{p-^p-^) < 00 such that for any 
u,v £ [—1, 1], any f E B, and any p £ M we have \g'p{u) — g'p{v)\ < C\u ~ v\. □ 

Proof of Theorem 14.21 The existence of solutions for 14.4|l has been proven in 
Lemma 14.31 to prove uniqueness we follow closely the argument in [Hj. Consider a 
solution F E T of 14.4|l . Since it solves H4.3|l almost everywhere, we get 

F'{u) = F'i-1) + J dwF'{w)n{p,F;w) (4.8) 

for all u in [—1, 1]. Moreover, taking into account that F is strictly increasing, we 
get from H4.3|l that 



expj 


[/ 


dwn{p, f;w) 


\ 


J dv exp 


{j 


f dwTlipJ; 


u;)| 



F' 

holds a.e., so that 

F{u)[l-F{u)] _ 



dv [1 - F{v) - p{v)] (4.9) 

1 



F'{u) F'{-1) 
for all u in [—1,1]. 

Let Fi, F2 G be two solutions of (|4.4II . If F{{-1) = F^i^l) an application of 
Gronwall inequality in (|4.8|) yields Fi = F2. We next assume F{{~1) < F2{—1) 
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and deduce a contradiction. Keep in mind that F/ > because Fi belongs to 
and recall H4.9|l . Let u :— inf{w G (—1, 1] : Fi{v) = F2{v)} which belongs to (—1, 1] 
because Fi(±l) = F2i±l) and F{{-1) < F^(-l). By definition of u, Fi{u) < F2{u) 
for any m e (— 1, u), i^i(u) = F2{u) and F[[u) > F2{u). By (|4.9|l . we also obtain 

Fi(zZ)[l-fi(ji)] F2(m)[1~F2(u)] 
F[{u) F^{u) 

or, equivalently, F[{u) < F2{u), which is a contradiction and concludes the proof 
of the first statement of Theorem 14.21 

We turn now to statement (i). Existence follows from identity (|4.3|l . which now 
holds for all points u in [—1,1] because p is continuous. Uniqueness follows from 
the uniqueness for the integro-differential formulation H4.4() . 

To prove (ii), let /?„ be a sequence converging to p in Al and denote by Fn = 
F{pn) the corresponding solution of (|4.4|l . By Lemma f4. 31 and Ascoli-Arzela the- 
orem, the sequence Fn is relatively compact in C^([— 1,1]). It remains to show 
uniqueness of its limit points. Consider a subsequence Uj and assume that Fn- 
converges to G in C^([— 1,1]). Since converges to p in and F„^. con- 
verges to G in C^([— 1,1]), by (|4.5|l ICp^.{Fnj) converges to JCp{G). In particu- 
lar, G ~ liiHj i^n^. — linij /Cp^ (F„^ ) — ^p{G) so that, by the uniqueness result, 
G — F{p). This shows that F{p) is the unique possible limit point of the sequence 
Fn, and concludes the proof of (ii). 

We are left to prove (iii). If p{t,u) e C^-°{[Q,T] x [-1,1]), we have from (i) 
and (ii) that F{t,u) = F{p(t, ■)){u) £ C°-^{[0,T] x [-1, 1]). We then just need to 
prove that F(t,u), as a function of t, is continuously differentiable. This will be 
accomplished by Lemma [4 . 41 b elow . □ 

In order to prove the differentiability of 1 1-^ F(t,u) :— F{p{t, •))(u) it is conve- 
nient to introduce the new variable 

Vit,u) :=log ^^^'^^ (i,^.)e[0,r]x[-l,l] (4.10) 

Note that ip e [(p_,(p+] where p± :— log[p±/(l — p±)] = log7-|- and u i-^ yy{t,u) 
is strictly increasing. We remark that, as discussed in 0], while the function F is 
analogous to a density, the variable ip can be interpreted as a thermodynamic force. 
The advantage of using p instead of F lies in the fact that, as a function of p, the 
functional Q is concave. This property plays a crucial role in the sequel. 

Let us fix a density profile p € (7^'° ([0,T] x [-1, 1]). By (i)-(ii) in Theorem Ol 
and elementary computations, we have that p G C°'^ ([0,T] x [—1, 1]) and it is the 
unique strictly increasing (w.r.t. u) solution of the problem 

A^P(h^+ ^-j-^ = pit,u) (t,u) e [0,T] X (-1,1) 

p{t,±l)^p± te[0,T] 

Note also that, by Lemma there exists a constant Ci = Ci(p_,p+) G (0, oo) 
such that 

7^ < V^(t , u) < Ci V (i , u) e [0, T] X [- 1 , 1] (4.12) 

Lemma 4.4. Let p E G^''^ {[0,T] x [—1,1]) and p — p{t,u) be the corresponding 
solution of !JlJ\ l. Then p G C^^^ qo,T] x [-1, 1]) and il){t,u) := dtp{t,u) is the 
unique classical solution of the linear boundary value problem 



(Vv?(t,M)) J (l + ev(*'"))- 



iP{t,u) ^ dtpit,u) (4.13) 
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for {t,u) £ [0,T] X (—1, 1) with the boundary condition 'ip{t,±l) — 0, t Cz [0,T]. 

Proof: Fix t G [0, T], for /i ^ such that t + h€ [0,T] let us introduce iphit, u) := 
[ifit+h, u)-ip{t, u)]/h. Note that, by (i)-(ii) in TheoremlO ■) € ([-1, 1]). 
By using 14.11|l . we get that iph solves 



VV'h(t, u) 



Wip{t,u)Wip{t + h,u)i (1 + e'^(*'")) (1 + e'^(*+'''")) h 
p{t + h,u) — p{t, u) 



for {t,u) e [0,T] X (-1, 1) with the boundary condition V'/ili, ±1) = 0, t G [0,T]. 

Multiplying the above equation by iphit, u) and integrating in du, after using the 
inequality x{e'^ — 1) > and an integration by parts, we get 



1 r^.i. u ..\\2 1-1 



, , /, . P{t + h,u) - p{t,u) 
du Iphit, u) 



h 



_i Vtp{t,u)Vtp{t + h,u) 

, , , ,7 ^ f\ { pit + h.u) - p{t,u)\ 



-1 

2 



where we used Schwartz inequality with e > 0. Recalling the Poincare inequality 
(with /(±1) = 0) 



du f'{u) 



2 



using H4.12|l and choosing e small enough we finally find 

limsup / du Iphit, u)]'^ < C2 / du [dtp{t,u)]'^ 

h^O J -I J-l 

for some constant C2 depending only on p^, 

Hence, by Sobolev embedding, the sequence iph{t,-) is relatively compact in 
C([— 1, 1]). By taking the limit /i — > in H4.14|l it is now easy to show any limit 
point is a weak solution of (|4.13() . By classical theory on the one-dimensional elliptic 
problems, see e.g. [22 IV, §2.1], there exists a unique weak solution of (|4.13(l which 
is in fact the classical solution since dtp{t, •) G C([— 1, 1]). This implies there exists 
a unique limit point ip{t, u) which is twice differentiable w.r.t. u. The continuity of 
t 1-^ ip{t, •) follows from the continues dependence (in the C^([— 1, 1]) topology) of 
the solution of (|4.13|l w.r.t. dtp{t, •) (in the C([-l, 1]) topology). □ 

The link between the boundary value problem H4.3|l and the variational problem 
(|2.11|) is established by the following theorem. 

Theorem 4.5. Let S be the functional on A4 defined in 1^2. Then S is bounded, 
convex and lower semi continuous on M. Moreover, for each p ^ M., we have that 
S{p) ~ Q(^p, F{p)^ where F{p) is the solution of \4-4}) - 

Proof: For each / G we have that G{-,f) is a convex lower semi continuous 
functional on Ai. Hence the functional S{-) defined in (|2.11|) . being the supremum 
of convex lower semi continuous functionals, is a convex lower semi continuous 
functional on Ai. Furthermore, by choosing f — p in H2.11|l we obtain that < 
Sa{p) < S{p). Finally, by using the concavity of x 1— > logs, Jensen's inequality, and 
/(±1) = p±, we get that Q{p,f) is bounded by some constant depending only on 
p_ and p+. 

In order to show the supremum in (|2.11l) is uniquely attained when / = F[p) 
solves (|4.4|) . it is convenient to make, as in Lemma the change of variables — 
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defined by ^(u) log - /(m)]}. Note that /(u) = e'^(")/[l + e'^(")]. 

We then need to show that the supremum of the functional 



du \ p{u) logp(u) + [1 - p{u)] log[l - p{u)] 



(4.15) 



+ [1 - p{u)]^{u) - log [1 + e^(")] + log 



for (/3 e := = {(^ e Ci([-l,l]j : ip{±l) = ip± , ip'iu) > 0} is uniquely 

attained when (f = (j){F{p)). We recall that F{p) denotes the solution of (|4.4|l . 

Since the real functions x ^ log a: and x ^ — log (l + e^) are strictly concave, 
for each p € M the functional g{p,-) is strictly concave on T. Moreover it is easy 
to show that g{p, •) is Gateaux differentiable on with derivative given by 



1 



1 + e¥'(") 



- p{u) 



By standard convex analysis, see e.g. ^2 I, Prop. 5.4], for any ip ^ if) € T we have 



g{p,^)<g{p,^) + (^ 



/sg{p, (fi) 



Sip 



By noticing that 5g{p, ip) / 5tf = Q \i Lp solves (j4.11|l a.e. we conclude the proof that 
the supremum on T of C/(p, •) is uniquely attained when tp — (f){F{p)). □ 



Remark 4.6. Given p E A4, let us consider a sequence pn E C^([— 1,1]) n Ai 
with p„(±l) ~ p±, bounded away from and 1, wtiicfi converges to p a.e. Then, by 
dominated convergence and (ii) in Theorem \4.^ we have S{pn) — g{pn, F{pn)^ — > 

g{p,Fip))^sip). 



5. The quasi potential 

In this section we show that the quasi potential for the one-dimensional bound- 
ary driven simple exclusion process, as defined by the variational problem H2.6|l . 
coincides with the functional S{p) defined in H2.11|l . In the proof we shall also 
construct an optimal path for the variational problem H2.5|l . 

Let us first recall the heuristic argument given in 4 . Taking into account the 
representation of the functional ItIttIp) given in Lemma 13.61 to the variational 
problem H2.6|l is associated the Hamilton-Jacobi equation 

where V denotes the derivative w.r.t. the macroscopic space coordinate u E [—1, 1]. 
We look for a solution in the form 

5V p f 
T- = log log ■ 



Sp ^l~p " 1 - / 
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and obtain a solution of (|5.1I) provided / solves the boundary value problem H4.3|l . 
namely / = F{p). On the other hand, by Theorem 14.51 we have 



ssip) _ sgipj) 



6p Sp 



, S5ipJ) 

f=F(p) 



SFjp) 

f=F(p) 



log log. 



l-p " 1 - F{p) 



since (j4.;^(l is the Euler-Lagrange equation for the variational problem (|2.1HI . We 
get therefore V = S since we have V{p) = S{p) = 0. 

Let TT*{t) = TT*{t, u) be the optimal path for the variational problem H2.6|l and de- 
fine p*{t) :— TT*{—t). By using a time reversal argument, in [Jj it is also shown that 
p*{t) solves the hydrodynamic equation associated to the adjoint process (whose 
generator is the adjoint of in L2{dp^)) which takes the form 

dtp*it) = -Up*it) +v(p*{t)[l - p*(t)]V^ ) (5.2) 

2 V op P=P'Wj 

We will not develop here a mathematical theory of the Hamilton-Jacobi equation 
(|5.1|) . We shall instead work directly with the variational problem H2.6|l . making 
explicit computations for smooth paths and using approximation arguments to 
prove that we have indeed V ^ S . Of course, the description of the optimal path 
will also play a crucial role. 

To identify the quasi potential V with the functional S we shall prove separately 
the lower bound V > S and the upper bound V < S. For this purpose we start with 
two lemmata, which connect S defined in H2.11|l to the Hamilton-Jacobi equation 
1)5. used for both inequalities. The bound V > S will then be proven by choosing 
the right test field H in ^T^. To prove V < S we shall exhibit a path Ti*{t) = 
iT*{t,u) which connects the stationary profile p to p in some time interval [0,T] 
and such that It{tt*\p) < S{p). As outlined above, this path ought to be the time 
reversal of the solution of the adjoint hydrodynamic equation (|5.2() with initial 
condition p. The adjoint hydrodynamic equation needs, however, infinite time to 
relax to the stationary profile p. We have therefore to follow the time reversed 
adjoint hydrodynamic equation in a time interval [0, Ti] to arrive at some profile 
p*{Ti), which is close to p if Ti is large, and then interpolate, in some interval 
[ri,Ti + Ts], between p*(Ti) and p. 

Recall that we are assuming p_ < and pick (5o > small enough for Sq < 
p_ < /5-i_ < 1 — Sq. For S e (0, Sq] and T > 0, wc introduce 

Ms := {p e f]) : p{±l) ^ p± , 5 < p{u) <l-5) (5.3) 



Dt,5 := {^eCi'2([0,T]x [-f,l]) : ^(t,±f)==p±, (5<^(<,u)<f-J} (5.4) 

Lemma 5.1. Let tt g Dt,s and denote by F{t,u) = F{-K{t, •)) {u) the solution of 
the boundary value problem {4--S^ with p replaced by 7r(t). Set 

, , 7r(t,u) F(t,u) , , 

F t, u) - log ^ \ > - log , y,/ , • 5.5 
1 — 7r(t, u) \ — t [t, u) 

Then, for each T > 0, 

S{n{T))-S{n{0))^ [ dt {dtn{t),r{t)) . (5.6) 
Jo 

Proof: Note that F{t, •) is strictly increasing for any t e [0, T] and F e C^'^ ([0, T] x 
[— f , 1]) by (iii) in Theorem l4.2l Moreover, since F{t, ±1) = yO±, we have dtF{t, ±1) = 
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0. By Theorem 14.51 dominated convergence, an explicit computation, and an inte- 
gration by parts, we get 

F{t) - 7r{t) AF{t) 



The lemma follows by noticing that the last term above vanishes by (|4.3|) . □ 

Lemma 5.2. Let p e Ais, denote by F{u) = F{p) (u) the solution of the boundary 
value problem i4.3)) , and set 

r{u) = log - log ■ 



1 - p{u) 1 - F{u) 

Then, 

(p(i-p),(vr)') + (Ap,r) = 0. (5.7) 

Proof: Note that F e Aig by Theorem 14.21 After an integration by parts and 
simple algebraic manipulations (|5.7|) is equivalent to 

We rewrite the first term on the left hand side as 



\ 'F{l-F)/ \ "F{l~F) 
which, by an integration by parts, is equal to 

/^r. \ / AF (l-2f)(Vf) 

\ ' F{l-F)/ \' 'F{l-F) [F(l-F)]2 
Hence, the left hand side of 1)5. 8|l is given by 

o-F^^\-/-^^l£— (n-F)' 
^ 'F(l-F)/ \[F(1-F)]2'^^ ' 

thanks to □ 

Note that, for smooth paths, Lemma l5.1l identifies. in the sense given by equation 
H5.6|l . r as the derivative of S. Lemma then states that this derivative satisfies 
the Hamilton-Jacobi equation (|5.1I) . 

5.1. Lower bound. We can now prove the first relation between the quasi poten- 
tial V and the functional S. 

Lemma 5.3. For each p G Ai we have V{p) > S{p). 

Proof: In view of the variational definition V, to prove the lemma we need to show 
that S{p) < /T(7r|p) for any T > and any path tt £ D(\0,T]; A4) which connects 
the stationary profile p to p in the time interval [0,T]: 7r(0) — p, 7r(T) — p. 

Fix such a path tt and let us assume first that tt g Dt,s- Denote by F(t) = 
F{Tr(t)) the solution of the elliptic problem (14.31) with 7r(t) in place of p. In view 
of the variational definition of It{tt\p) given in H2.5|l . to prove that S{p) < /T(7r|p) 



21 



it is enough to exhibit some function H E Co'^([0,T] x [-1,1]) for which S{p) < 
Jt,h,p{t^)- We claim that F given in (|5.5() fulfills these conditions. 

We have that T e Co'^([0,r] x [-1, 1]) because: tt e Dt.s by hypothesis, F e 
C^'2([0,r] X [-1,1]) by (iii) in Theorem Ol r(t,±l) = since 7r(t, •) and F{t,-) 
satisfy the same boundary conditions. Recalling (|2.4|l we get, after integration by 
parts, 

Jt.vA^) - / dt{dtTr{t),T{t)) 
Jo 

- l£dt [(r(i),A^(t)) + (4t)[i-^(t)],[vr(i)]')] . 

By Lemmata 15. II and 15. 21 we then have JT,r,p{^) — S{p)- 

Up to this point we have shown that S{p) < IxiT^lp) for smooth paths tt bounded 
away from and 1. In order to obtain this result for general paths, we just have 
to recall the approximations performed in the proof of the lower bound of the large 
deviation principle. Fix a path tt with finite rate function: /^(Trlp) < oo. In Sec- 
tion |^| we proved that there exists a sequence {7r„, n > 1} of smooth paths such 
that 7r„ converges to tt and /T(7r„|p) converges to It{t^\p)- Let 7f„ be defined by 
(1 — n^^)7r„+?i~^p. Since 7r„ converges to tt, 7r„ converges to tt. By lower semi con- 
tinuity of the rate function, /^(Trlp) < liminf„^oo /^(Trnlp). On the other hand, by 
convexity, /T(7f„|p) < {1 - n-^)lT{Trn\p) + n'^ It{p\p) = (1 - n"^)/T(7r„|p) so that 
limsup„^oQ /T(7r„|p) < /T(7r|p). Since nn belongs to Dt,s for some S = Sn > 0, each 
path TT with finite rate function can be approximated by a sequence 7r„ in Dt,s„ , for 
some set of strictly positive parameters Sn, and such that /T(7r|p) = lim„ /T(7^n|p)■ 
Therefore, by the result on smooth paths and the lower semi continuity of S, we 
get 

iTiMp) = lim/T(^„|p) > hminf 5(^„(T)) > S{tt{T)) 

n n 

which concludes the proof of the lemma. □ 



5.2. Upper bound. The following lemma explains which is the right candidate 
for the optimal path for the variational problem (|2.6f) . 

Lemma 5.4. Fix 6 G (0,(5o], a profile a G Ms, o,nd a path tt G Dt,5 with finite 
rate function, IxiTrla) < oo. Denote by F{t,u) — F{Tr{t, •)) (u) the solution of the 
boundary value problem with p replaced by 7r(t). Then there exists a function 

K G 7ii(7r) such that tt is the weak solution of 

' 5t7r= -iA7r + v(^(l-7r)V[log^i^ + if]) (t, u) G [0, T] x (-1, 1) 

< ^(t,±i) = p± te[o,r] (5.9) 

7r(0, u) = a(u) u G [—1, 1] 

Moreover, 

/T(^|a) = S(it(T)) - S{ot) + i ^\t {-Km - ^(t)], \^K{t)f) (5.10) 

The optimal path for the variational problem (|2.6|l will be obtained by taking a 
path TT* for which the last term on the right hand side of the identity H5.1()|l (which 
is positive) vanishes, namely for a path tt* which satisfies H5.9(l with if = 0. Then 
p*(t) = •7r*(— t) wiU be a solution of H5.2|) . 

Proof: Denote by Fl the function in Hi{tt) introduced in Lemma 13.61 let F as 
defined in H5.5() . and set A' := F — H. Note that K belongs to 7ii(7r) because: 
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TT e Dt,5 by hypothesis, F G C^'^{[0,T] x [-1, 1]) by Theorem lO and T{t,±l) = 
0. Then 1)5.9(1 foUows easily from 1(8. To prove the identity ((5. 10(1 . replace in 
1(5.6(1 dtTr{t) by the right hand side of the differential equation in ((5.9() . After an 
integration by parts we obtain 

SiniT))~Sia) = /"^d Ji(r(t), A^(t)) + (7r(i)[l - ^(t)], [Vr(i)F) 







where we used Lemma [5 .21 Recalling K = T — H, we thus obtain 



1 '■^ 



2./0 







which concludes the proof of the lemma in view of 13.4(1 . □ 

We write more explicitly the adjoint hydrodynamic equation ((5.2(1 . In the present 
paper, we shall use it only to describe a particular path which will be shown to be 
the optimal one. For p £ M, consider the non local differential equation 

dtp* = ^Ap* - V(p*(l - p*)V log ^) (i, u) e (0, oo) X hi, 1] 



F{t,u) = Fip*{t,-)){u) (i,M)e (0,cx3)x [-1,1] 

p*(t,±l) = P± te{0,(x) 
P*i0,u) = p{u) ue[^i,i] 



(5.11) 



where we recall that F{t, u) = F{p*{t, •)) (it) means that F{t, u) has to be obtained 
from p*{t,u) by solving 1(4.4(1 with p(u) replaced by p*{t,u). Since Vlog[i^/(l — 
F)] > 0, in 1(5.11(1 there is a positive drift to the right. Let us describe how it is 
possible to construct the solution of ((5.11(1 . 

Lemma 5.5. For p (E A4 let $(i) be the solution of the heat equation \2.1'2\) and 

define p* ^ p *{t,u) by ^J^) . T/ien p* G C^-^ ((0, oo) x [-1, 1]) n C([0, oo); 7W) and 
solves 1^5. Moreover, if S < p{u) <1 — S a.e. for some 6 > 0, there exists S' = 
(5'(p_,p+,(5) e (0,1), for which 6' < p*{t,u) < 1-5' for any {t,u) e (0, oo) x [-1, 1] . 

Proof: Let F{u) = F{p)(u), then, by Theorem lO F e Ci([-l,l]) and, by 
Lemma 14.31 there is a constant C G (0, oo) depending only on p_, p^ such that 
(7-1 < F'{u) < C for any u e [-1,1]. Since ^{t,u) solves ((2.12(1 . there exists 
Ci = Ci(p_,p+) e (0,oo) such that Cf^ < (V$)(i,it) < Ci for any {t,u) e 
[0,oo) X [-1,1]. Moreover, $(t,±l) = p± so that A$(i,±l) = 2dt^{t,±l) = 0. 
Hence, p* defined by 1(2.13(1 satisfies the boundary condition p*{t, ±1) — $*(t, ±1) = 
p±. Furthermore p* e C^'^((0,oo) x [—1,1]) . 

For the reader's convenience, we reproduce below from T, Appendix B] the proof 
that p*{t,u), as defined in ((2.131) . solves the differential equation in ((5.11(1 . From 
(j^T^ we get that 

' 'i' - 1 + (1 - 2$)^=^ - $(1 - $)^ 

(V$)' (V$)' 
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recalling (|2.12l) , by a somehow tedious computation of the partial derivatives which 
we omit, we get 



2 



(V$) J V$(l-$) 

from which, by using again H2.13|l . we see that p* satisfies the differential equation 
in ifTTTIl . 

To conclude the proof of the lemma, notice that p* is the solution of 
dtp* ^\/\p* ~^{p*{l-p*)VH} , 
p*(t,±l)-p± , 
P*(0,-)=P(-), 

for some function H in C^'^ ([0, oo) x [—1, 1]) for which Vi? is uniformly bounded. 
Though H does not vanish at the boundary, we may use a weakly asymmetric 
boundary driven exclusion process to prove the existence of a weak solution A(t, u), 
in the sense of Subsection l3.4l which takes values in the interval [0, 1]. Since VH is 
bounded, the usual method gives uniqueness so that \ — p* and < p* < 1. 
In particular p* e C([0, oo); 7W) . 

Assume now that 5 < p < 1 — (5 for some 6 > 0. Fix t > and assume that 
p*{t, •) has a local maximum at —1 < uq < 1. Since p* is a smooth solution of 
H5.11(l . a simple computation gives that at (t, uq) 

{dtp ) = -Ap F^l-F)^ +F-1) 

because (Vp*)(t,uo) = and Alog{F/l-F} = {VF)^{p*+F-l)/F^{l~Ff. Since 
Mo is a local maximum, Ap* < 0. On the other hand, assume that p*{t, uq) > 1 — p_, 
in this case, since p_ < _F, p* + F — 1 > so that dtp* < 0. In the same 
way we can conclude that {dtp*){t,ui) > if ui is a minimum of p*{t,-) and 
p*(t, Ml) < 1 — p+. These two estimates show that min{5, 1 — p_|_, p_} < p*{t, u) < 
max{l — 5, 1 — p_, p-|_}, which concludes the proof of the lemma. □ 

We now prove that the solution of (|5.11(l . as constructed in Lemma l575l converges, 
as t ^ oo, to p uniformly with respect to the initial datum p. We use below the 
usual notation ||/||oo := supug[_i,i] |/(u)|. 

Lemma 5.6. Given p G M, let p*{t) — p*{t,u) be the solution 1^,5. Then, 

lim sup \\p*{t) - p\\ = . 

Proof: Let us represent the solution $(i) of H2.12|l in the form u) = p{u) + 
^{t,u). Then *(i) = P°^'(0) where P° is the semigroup generated by (1/2)A°, 
with A" the Dirichlet Laplacian on [—1, 1]. Since ^'(0) — F{p) — p and since the 
solution F{p) of (|4.4() as well as p are contained in the interval [p_,p+], we have 
that ||^'(0)||oo < |p+ — P-l < 1- Therefore, by standard heat kernel estimates, 

lim sup |||*(i)||oo + ||V*(t)||co + ||A*(<)||eo| = 

the lemma follows recalling that, by Lemma [5. 51 p*(t) is given by (|2.1,'^|l . □ 

Lemma |5.6I shows that we may join a profile p in to a neighborhood of the 
stationary profile by using the equation 15.11|l for a time interval [0, Ti] which at 
the same time regularizes the profile. On the other hand, from Lemma 15.41 we shall 
deduce that this path pays S{p) — 5(p*(Ti)). It thus remains to connect p*(T'i), 
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which is a smooth profile close to the stationary profile p for large Ti , to p. In the 
next lemma we show this can be done by paying only a small price. We denote by 
II • 1 1 2 the norm in L2 ( 1 , 1] , du) . 

Lemma 5.7. Let a G Msa a smooth profile such that \\a — p||oo < 

Then there exists a smooth path Tvit), t G [0, 1] with 5q/2 < tt < 1 — <5o/2, namely 

n £ with 7r(0) = p, tt{1) = a and a constant C = C{So) € (0, 00) such that 

h{n\p)<C\\a-p\\l 

In particular V{a) < C\\a — p\\2- 

We remark that by using the "straight path" 7r(i) = p(l — t)+at one would get 
a bound in terms of the Hi norm of a — p. Below, by choosing a more clever path, 
we get instead a bound only in term of the L2 norm. 

Proof: Let {ek,Xk), A: > 1 be the spectral basis for — (1/2)A°, where A'^ is the 
Dirichlet Laplacian on [—1, 1], namely {efc}fc>i is an complete orthonormal system 
in L2{[— 1,1], du) and — (l/2)A°efc = AfeCfc. Explicitly we have ek{u) = cos(fc7ru/2) 
and Afe = We claim that the path n{t) = n{t,u), {t,u) e [0,1] x [-1,1] 

given by 

7^(0 =P + £^>^— {(«-P,efe)efe (5.12) 
fc=i 

fulfills the conditions stated in the lemma. 

It is immediate to check that 7r(0) = p, ^"(1) = a and 7r(t, ±1) = p±. Further- 
more, by the smoothness assumption on a, we get that tt € C^'^[[Q, 1] x [—1, 1]). 
In order to show that (5o/2 < tt < 1 — 5q/2, let us write ^{t) = p + q{—t), then 
q{t) —q{t,u), {t,u) e [—1,0] X [-1,1] solves 

dtq{t) = iAg(i) + g 

q{t,±l) = 
q{-l,u) = a{u) - p{u) 

where g = g{u) is given by 

k=l 

Let us denote by ||.g||ffi := ||.9'||2 the Hi norm in [—1,1]; a straightforward 
computation shows 

00 \ 2 s °° 

IMIl,^ = ^2A.(-3^) (a-p-,e.)^<^^(„_p-e,)2 

k=l ^ ^ k=l 

2 '<5o\2 1 ,2 



where we used that, for A > 0, we have e'^ — 1 > A^/2. 

Let P^ = exp{iA''/2} be the heat semigroup on [—1, 1]; since ||g||oo < V2||5||hi 
we have 



9(0 II 00= sup P°_^i(a-p)+ f dsP^_^g 
tef-1,01 J-i 



sup 

te[-i,o] te[-i,o] 



16 TT 16 



so that 7T g £'i,5(,/2- 

We shall estimate /i(7t|p) by using the representation given in Lemma 13.61 To 
this end, let us define h = h{t,u) e C([0, 1] x [-1, 1]) by h := -dtir + {1/2)At: and 
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let H ~ H{t, u) be the solution of 

V {tt[1 - tt]V H) =/i 
H{t,±l) = 

so that TT solves with H as above which belongs to Hi{n). 

Let us denote by || • the usual negative Sobolev norm in [—1, 1], namely 

I|/^IIL.:= sup j^j^-f:^{Ke,r 
By using that 7r[l — tt] > ((5o/2)^ a simple computations shows 

\t {nm mi (vi?w)') < ^ fdt \\h{t)\\i_^ 

By using the explicit expression for tt we get 

fc=i 

hence, by a direct computation, 



oo 



where we used that for A > Ai we have e > 2. We thus get 



fe=i 

A 



2 /■ 1 8 ^ 

/i(^r|p)< ;^ / dt||/^(t)||Li < 72E' 



2__2_|l -||2 



a-p,ek) ^j2\\^~P\\2 



k=i 

which concludes the proof of the lemma. □ 

We can now prove the upper bound for the quasi potential and conclude the 
proof of Theorem 12. II 

Lemma 5.8. For each p ^ Ai, we have V{p) < S{p). 

Proof: Fix < e < (5o/(32), p e M and let p*{t,u) be the solution of (|5.11(l with 
initial condition p. By Lemma l5^ there exists Ti = Ti(e) such that ||p*(0~Plloo < ^ 
for any t >Ti. Let a := p*{Ti) and let tt be the path which connects p to a in the 
interval [0, 1] constructed in Lemma 15.71 

Let T := Ti + 1 and ■n*{t), t e [0, T] the path 

_ r m forO<t<l 
W - I p*(T-i) for 1 <t <T ^^•^'^^ 

By Remark 14 .61 given p e as above, we can find a sequence n > 1} with 
p„ G M-S„ for some (5„ > converging to p in Al and such that S{pn) converges to 
S{p). Let us denote by p"'* the solution of (|5.11() with initial condition p„ and set 

^ > \ p''^*{T -t) for 1 < i < T ' 

where 7r"'*(i) is the path joining p to a„ := p"'*(Ti) in the time interval [0,1] 
constructed in Lemma [5.71 We claim that the path tt"'* defined above converges 
in D([0,T], Al) to TT*, as defined in (|5.13|) . Before proving this claim, we conclude 
the proof of the lemma. 
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By the lower semi continuity of the functional It{-\p) on D{\0,T], Ai) we have 
It{tt*\p) < liminf /T(7r"'*|p) (5.15) 

n 

On the other hand, by definition of the rate function and its invariance with respect 
to time shifts we get 

/T(7r"'*|p) =/i(7r"'*|p)+/T,(p"'*(Ti-.)|p"'*m)) (5.16) 

By Theorem Ol F„ := F(p„) converges to F = F{p) in C^{[-l,l]) so that 
$„(t), the solution of l|2.9|l with initial condition Fn, converges to ^{t) in ([—1, 1]) 
for any t > 0. Hence, by ^J^, /9"'*(Ti) converges to p*{Ti) in C([-l, 1]). Recall- 
ing that ||p*(Ti) — p\\oc < £ < i5o/(32), we can find A^o = -^o(<5o) such that for any 
n > A^o we have ||p"'*(Ti) — p||oo < ^ < <5o/(16). We can thus apply Lemma [5.71 
and get, for n> Nq 

/i(7r"^*|p) <C|1p"^*(Ti)-p||2 (5.17) 

for some constant C — C{So). 

By Lemma [5.51 p"'*(Ti — t), t E [0, Ti] is smooth and bounded away from and 
1, namely it belongs to -Dti,i5„ for some (5„ > 0. We can thus apply Lemma and 
conclude, as p"'*(Ti - t) solves (jSSl) with K ^ 0, 

/t,(p"^*(Ti - ■)\p-'*{T,)) = S{pn) - ^(p"'*(ri)) < S{p^) (5.18) 

From equations (|5.15|) - H5.18(I we now get 

lT{7r*\p) < hminf [^(p„) + C|lp"'*(ri)-p|l^] 

n 

= Sip) + C\\p*{T,)-p\\l<S{p) + 2Ce^ 

and we are done by the arbitrariness of e. 

We are left to prove that tt"'* — s- tt* in D([0, T],M). We show that tt"'* converges 
toTT* inC{[0,T];M). Pickei e (0, Ti]; since p"^*(t) converges to p* (i) inC([-l,l]) 
uniformly for t € [ei, Ti] we conclude easily that tt"^* converges to tt* in C([l, T — 
El] X [—1,1]). We recall that, by Lemma lOI VFn{t) and VF(t) are uniformly 
bounded. Moreover, 7r"^*(T - t) and tt*{T -t), t e [T -Ti,T] are weak solutions 
of (|5.11|) : for each G G C([-l, 1]) we thus get 

lim limsup sup |(7r"^*(i), G) - (7r*(t), G)| = 

'^iJ-O n te[T-ei,T] 

this concludes the proof that p"'* converges to p* in C{\l,T\]M). Since p*'"(ri) 
converges to p*'"(Ti) in G^([— 1, 1]) it is easy to show that tt"'* converges to tt* in 
C([0, 1] X [-1,1]). Hence tt"'* converges to tt* in C([0, r];7W) □ 



Appendix A. A lower bound on the quasi potential {d > 1) 

In this Appendix we prove a lower bound for the quasi potential in the d dimen- 
sional boundary driven simple exclusion process. For d = 1 this bound has been 
derived from l|TTT|l in 001 • 

Let A C M'' be a smooth bounded open set and define Ajv := n A^A. Let also 
7(m) be a smooth function defined in a neighborhood of 9A. The d-dimensional 
boundary driven symmetric exclusion process is then the process on the state space 
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Sat := {0, 1}^" with generator 
2 



{x,y}CAN 
\x-y\ = l 

+ ^ E ('?(^) + [1 - )) - fiv)] 

ajSAiv.y^Ajv 
2;-y| = l 

where tr^'^ and cr^ have been defined in Section |21 

The hydrodynamic equation is given by the heat equation in A, namely 

dtp = iAp ue A 
p{t, u) = a{u) u e i9A 
p{0,u) = po{u) 

where a{u) — ^{u)/[l + j{u)]. We shall denote by p = p{u), u £ A the unique 
stationary solution of the hydrodynamic equation. 

By the same arguments as the ones given in Section O it is possible to prove the 
dynamical large deviation principle for the empirical measure. The rate function is 
still given by the variational formula H2.5|l . but we now have 

Jt,hA^) {tt{T),H{T))-{p,H{0))- I dt (^7r{t),dtH{t) + ^AH{t)) 

1 , .2. 1 '•^ 



dt{x{TT{t)),{VH{t)) ) + - / dt d<j{u) a{u)dnH{t,u) 

^ Jo ^ Jo JdA 

where dnH{t,u) is the normal derivative of H{t,u) {fi being the outward normal 
to A) and a{u) is the surface measure on dA. 

Let us define the quasi potential V{p) as in H2.6|l and set 



So{p) := / du 



piu)log^ + [l- p{u)]\oi 



p{u) l~p{u)_ 
Theorem A.l. For each p £ A4 we have V{p) > So{p). 

Proof: We shall prove that It{tt\p) > So{p) for any 7r(-) such that 7r(0) = p 
and 7r(T) = p. Let us assume first that tt G C^''^{[0,T] x A), 7r(i, u) = a{u) for 
{t, u) G [0, T] X dA, and vr is bounded away from and 1. Given such tt we use the 
variational characterization of It and chose 

X , '!r(t,u) , p(u) 
H{t,u)= log ^ ^'./. -log- ' 



1 — 7r(<, u) 1 — p{u) 

Note that Hit, m) = for [t, u) e [0, T] x dA since tt and p satisfy the same boundary 
condition. By dominated convergence and an explicit computation we get 

5o(^(r)) - SoiT^m = Cdt j^So{n{t)) = Cdt {dtTT{t),H{t)) 



Recalling that Jt.h,p{'^) has been defined above, a simple computation shows 
Jt.hA^^) - SoiAT)) + \ I dt {\/H{t),\/n{t) - n{t)[l - TT{t)]\/ H{t)) 



2 _ „ 

= 'oi<T)) + l£dt{^^^i.it)-pf 

since the second term above is positive we conclude the proof of the lemma for 
smooth paths. To get the general result it is enough to repeat the approximation 
used in Lemma l5.3l □ 
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